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The spherical to gamma-unstable nuclei shape- phase transition in odd-A nuclei is investigated 
by using the Dual algebraic structures and the affine 517(1,1) Lie Algebra within the framework of 
the interacting boson - fermion model. The new algebraic solution for A-odd nuclei is introduced. 

In this model, Single j = 1/2 and 3/2 fermions are coupled with an even-even boson core. Energy 
spectra, quadrupole electromagnetic transitions and an expectation value of the d-boson number 
operator are presented. Experimental evidence for the U BF (5) — O bf ( 6) transition in odd -A Ba 
and Rh isotopes is presented. The low-states energy spectra and B(E2) values for these nuclei have 
been also calculated and compared with the experimental data. 
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I. INTRODUCTION 

Quantum Phase Transitions ( QPTs ) are sudden changes in the structure of a physical system. Nuclear physics has 
important contributions to make to their study because nuclei display a variety of phases in systems ranging from few 
to many particles [lj.The signs of QPTs in nuclear physics are changes in mass and radius of nuclei and quantities 
such as level crossing and electromagnetic transition rates when the number of protons or neutrons is modified. Phase 
transition happens in both even-even and odd — A nuclei. Phase transitions investigations have been mostly performed 
on even-even systems 0 within the framework of the interacting boson model {IBM) ii ,which describe nuclei in 
terms of correlated pairs of nucleons with L = 0,2 treated as bosons ( s 7 dbosons ). The IBM Hamiltonian has exact 
solutions in three dynamical symmetry limits (17(5),517(3) , and 0(6)). These situations correspond to the spherical, 
axially deformed, and gamma-unstable ground state shapes, respectively. The transition between the 17(5)-517(3) 
limits is a first-order shape-phase transition while a second-order shape-phase transition occurs between the 17(5) 
and 0(6) limits [H-[E] . During a transition from one limit to another, meet the points in which potential has flat 
behavior. These points are called critical point.Recently Iacliello introduced the so-called critical point symmetries in 
the framework of the collective model for even-even nuclei. The critical point from spherical to 7 — unstable shapes, 
called E{ 5) [l,|l|, the critical point from spherical to axially deformed shapes, called X (5) [3, 6j, and the critical point 
from axially deformed shapes to triaxial shapes, called Y(5) @,|7]. Phase transitions is also investigate in odd-A nuclei 
within the framework of interacting boson - fermion model ( IBFM )[§[, which describe nuclei in terms of correlated 
pairs, with L = 0, 2 (s , d bosons), and unpaired particles of angular momentum j {jfermions ).Studies of QPTs in 
odd-even nuclei were implicitly initiated years ago by Scholten and Blasi Q. Several explicit studies have recently 
been made by Alonso et al. floi - ITil ] and by Boyukata et al. 0 , who also have suggested a simple form of the IBFM 
Hamiltonian, particularly well-suited to study QPTs in odd-even nuclei because of their supersymmetric properties. 
Similar to even-even nuclei, also exists critical point for odd-even and even-odd nuclei but in this case, critical points 
show with E{5/'}2j 2 j + 1) and X(5/JT 2 j + 1) that j is the angular momentum of single nucleon. Iacliello [lj, [lj| 
has been the case of a j(3/2) fermion coupled to a boson core that undergoes a transition from spherical to 7 — unstable 
shapes. At the critical point, an elegant analytic solution, called E{ 5/4), has been obtained starting from the Bohr 
Hamiltonian fl4j |. 

In this study, we investigate the transition U BF { 5) — O bf { 6) in odd — A nuclei. The new algebraic solution for 
A-odd nuclei is introduced. For this transition only the boson core experience the transition and fermion with j = 1/2 
and 3/2 coupled to boson core. We evaluate exact solutions for eigenstate and energy eigenvalues for transitional 
region in the IBFM by using the Dual algebraic structure for the two level pairing model that based on Richardson 

- Gaudin method and changing the control parameter that based on affine 517(1,1) lie algebra. In order to the 
investigation of phase transition, we calculate observables such as level crossing, expectation values of the d-boson 
number operator, ground - state energy and its first derivative . The low-lying states of 56 Ba and 45 ~~ Rh 
isotopes have been studied within suggested model. The results of calculations for these nuclei will present for energy 
levels and transitions probabilities, two neutron separation energies and will compare with the corresponding the 
experimental data. 

This paper is organized as follows: Section 2 briefly summarizes theoretical aspects of transitional Hamiltonian and 
affine 517(1,1) algebraic technique. Sections 3 and 4 include the results and Experimental evidence and sect. 5 is 
devoted to the summary and some conclusions. 


II. THEORETICAL FRAMEWORK 

The 517(1,1) Algebra has been explained in detail in Refs p~6j - [l8| .The 517(1,1) algebra is produced by S u , v = 0 
and ± , which satisfies the following commutation relations 

[5°,5 ± ] = ±5 ± , [5+5“] = -25° (1) 

The quadratic Casimir operator of 517(1,1) can be written as 

C 2 = 5° (5° - 1) - 5+5~ (2) 

The basis states of an irreducible representation ( irrep) 517(1,1}_,|are determined by a single number k , k where 
can be any positive number and fi = k,k + 1,... Therefore fl71. Iisj|. 


02(517(1,1 ))M = Kk - 1)| kfi) 


S°\k/j,) = n\kfj) 


(3) 


3 


In IBM , the generators of SU d ( 1,1) generated by the d-boson pairing algebra 

S + (d) = i(d+.d+) ,S~(d) = ^(d.d), S\d) = -J2{ d i d v + d v d i) = \Vn d + 5) (4) 

V 

Similarly, s- boson pairing algebra forms another SU S (1, 1) algebra generated by 

5' + (s) = ^ s+2 » S ~( s ) = ^s 2 , 5°(s) = i(s+s + ss+) = i(2n g + l) (5) 

SU sd (l, l)is the s and d boson pairing algebras generated by 

S + (sd) = ~(d + .d + ± s + ) , S~(sd) = —(d.d ± s 2 ) , 5°(sd) =+ dj,d+) +-(s + s + ss + ) (6) 

V 

Because of duality relationships El, [2(j, It is known that the base of U( 5) D SO( 5) and 50(6) D SO( 5) are 
simultaneously the basis of SU d ( 1,1) D U( 1) and SU sd ( 1,1) D 17(1), respectively. By use of duality relations [l7lfl9l], 
the Casimir operators of 50(5) and 50(6) can also be expressed in terms of the Casimir operators of SU d ( 1,1) and 
SU sd ( 1,1), respectively 


0 2 (5[/ d (l,l)) = ^ + i(7 2 (50(5)) 


(7) 


0 2 (5t/ sd (l,l)) = ^ + i(7 2 (50(6)) 


The infinite dimensional 517(1,1) algebra that is generated by use of EM 

, S° n =c 2n S 0 (s) + c 2n S 0 (d) 


Sn = c s n+ 1 5 ± (s) + c d n+1 S ± (d) 


( 8 ) 

(9) 


Where c s and Cd are real parameters and n can be 0, ±1, ±2,. These generators satisfy the commutation relations 


[51, St] = ±St+n , [S+, 5-] = -25^ + „ +1 


( 10 ) 


Then,5^, fi = 0, +, —; m = ±1, ±2, ...generate an affine Lie algebra517(l, l)without central extension. 

In odd A nuclei the Bose — Fermi symmetries are associated with each of the dynamic symmetries of IBM-lQ. 
So, the boson algebraic structure will be always taken to be U B ( 6 ) , while the fermion algebraic structure will depend 
on the values of the angular momenta, j, taken into consideration [§]. First we considered the case that a system of 
N bosons (with L=0, 2) coupled to a fermion with angular momentum j=l/2. The Lattice of algebras in this case is 
shown in Fig. 1.In Figs 1, 2 , the chain 1 show the state that bosons have U B (5) dynamical symmetric while bosons 
in chain 2 have O b ( 6 ) dynamical symmetric.By employing the generators of Algebra 517(1,1) and Casimir operators 
of subalgebras , the following Hamiltonian for transitional region between U BF (5) — O bf (6) limits is prepared 


H = gS+So + aSi + pC 2 {SO B (5)) + SC 2 {SO B (3)) + 1 C 2 (spin BF (3)) (11) 

Following this, we considered the state that the odd nucleon being in a j=3/2 shell. The Lattice of algebras in this 
case is also shown in Fig.2 .The following Hamiltonian for state that odd nucleon being in a j=3/2 shell for transitional 
region between U BF (5) — O bf ( 6 ) limits is prepared 


H = gS+So + aS ? + pC 2 {s P in BF (5)) + 'yC 2 (spin BF (3)) (12) 

Eqs.(ll) and (12) are the suggested Hamiltonians for boson - fermion systems with j = 1/2, 3/2, respectively and 
a , j3 ,6 ,7 are real parameters. By considering Eqs.(2),(7),(8), It can be shown that Hamiltonians (11) and (12) 
are equivalent with O bf (6) Hamiltonian when c s = c d and with U BF (5) Hamiltonian if c s = 0 and c d ^ 0. Thus, 
As mentioned because that only the boson core experience the transition and fermion coupled to boson core , the 
c s 7 ^ c d ^ 0 situation just corresponds to U BF ( 5) «-»■ O bf (6) transitional region. In our calculation, we take c d {= 1) 
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U B (6) 0 u r (i) 



Spin BF ( 2) 


FIG. 1. The Lattice of algebras in the case that a system of N bosons (with L = 0,2) coupled to a fermion with angular 
momentum j = 1/2. 


U b (6') 0 U f (4) 



FIG. 2. The Lattice of algebras in the case that a system of N bosons (with L — 0, 2) coupled to a fermion with angular 
momentum j = 3/2. 


constant value and c s change between 0 and Cd■ For evaluating the eigenvalues of Hamiltonians (11) and (12) the 
eigenstates are considered as 


| k; v s vni\LM) = 


E 


Q+ C+ C+ Q+. 


ln 1 n 2 ....n k X 1 X 2 X 3 ...X k d ni b n 2 d n3 - 


rii(zZ 


(13) 


Eigenstates of Hamiltonians (11) and (12) can obtain with using the Fourier-Laurent expansion of eigenstates and 
SU( 1,1) generators in terms of c-unknown number parameters Xi with i = 1, 2,..., k. It means, one can consider the 
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eigenstates as 0,0 

\k; v s vn^LM) = JVS+5+5+...5+ \lw) BF 
Where N is the normalization factor and 

S+ = . E S+{s) + 1 E S+(d) 


1 - C 2 s Xi 


1 ^ c d X i 


The c-numbers Xi are determined through the following set of equations: 

a _ gc 2 (u s + \) + gc 2 d (y d + f) _ ^ 2 g 


Xi 1 - C 2 Xi 


1 - c d Xi *rf.Xi- Xj 
a 


(14) 

(15) 

(16) 


With Clebsch — Gordan(CG) coefficient, we can calculate lowest weight state, \lw ) BF , in terms of boson and fermion 
part.For the j=l /2 case we have: 


1 5 

Mm±i = I N ,k d = 2 ( u d + 2 ).W = + ^)>Ms = + A),£,m± A) (17) 


^(n,+J 


1 , 1 , 
2 O'* + 2 ^ 


2 Os + 2 ^ 


\ R P / T i Til Q O / T -p 111’ Q O 

w ) = ± \/ /or , in M m ix+ + y /o r , ^ IM m+ ix- 


(2 L + 1) 

The lowest weight state for the j=3/2 case is calculated as: 


1 


1 


(2L + 1) 1 /m+ 2' 


1 , 1 , 


> m ±! = 1^,0 = gO d + 2^ ,Md = ^Orfd- 2)Os = 2Os + 2)-/^s = gO* + 2)’ i > m± 2^ 

= §."»j = !) + C^_ 4 , 4 |Mf ili = |,mj = I) 


y) BF = C 
' 1 


m,m— 
3 


3 3 | IW) 
2 ’ 2 


+ (T’ t , \lw) B + x |j = §, mj = -h) + <7 


J,L, 
m,m 

• 7 > L >f 17.. AS 


2 5 2 
3 


m , m+ 3 j _3lM m+ |b'= 2 ,mj = - 2 > 


The C r OrO,m,- symbols represent Clebsch-Gordan coefficients, where 

EM bf = a°M bf 

A° =cO(^ + i)i+c^(^ + |)i 
The eigenvalues of Hamiltonians (11), (12) can then be expressed; 

£ (fc) = O fe) + aA° + 0i/ d (i/ d + 3) + SL(L + 1) + 7 J(J + 1) 

£ (fc > = /i (fc) + «A° + ^(i/i(i/i + 3) + ^2 + 1)) + 7 J( J + 1)) 


(18) 

(19) 

( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 


h(k) = E 

i= 1 


a 


Xi 


(25) 


The quantum number (/c) is related to the total boson number N by 

N = 2 fc -f i^ s + 

In order to obtain the numerical results for energy spectra ( E of considered nuclei, a set of non-linear Bethe — 
Ansatz equations (BAE) with k- unknowns for k-pair excitations must be solved 00 also constants of Hamiltonian 
with least square fitting processes to experimental data is obtained. To this aim, we have changed variables as 

G = — <l,g = l,yi = c 2 d Xi 

Cd 
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So, the new form of Eq.(16) would be 

“=^L±il + ^±i 2 -£_L_ (26) 

Vi 1 - C-yi I-yi ^ Hi - Vj 

To calculate the roots of Bethe — Ansatz equations (BAE) with specified values of v s and Vd, we have solved Eq. (24) 
with definite values of C and a 0 .Then, we carry out this procedure with different values of C and a to give energy 
spectra with minimum variation in compare to experimental values filll ; 

wtot T-r. 
l,tot 

( N to t the number of energy levels where included in the fitting processes). The method for optimizing the set of 
parameters in the Hamiltonian (/3,y ,6) includes carrying out a least-square fit (LSF) of the excitation energies of 
selected states 0 


III. RESULTS: 

This section presented the calculated phase transition observables such as level crossing, ground - state energy and 
the derivative of the energy, expectation values of the d-boson number operator and energy differences. 


A. energy spectrum and level crossing 

To display how the energy levels change as a function of the control parameter C and the total number of bosons 
N, the lowest energy levels as a function of C for N = 10, 20 bosons are shown in Fig.3 and Fig.4, where in Fig.3 
other fixed parameters are a = 1000, /3 = —57, <5 = 41, 7 = 36 and Fig.4 obtained with a = 1000, /3 = 6.5, 7 = —22. 
Figs show how the energy levels as a function of the control parameter C evolve from one dynamical symmetry limit 
to the other . It can be seen from Figs that numerous level crossings occur; especially in the region around C > 0.7. 
The crossings are due to the fact that i/d, 0(5) quantum number called seniority, is preserved along the whole path 
between 0(6) and 17(5) (22j. With increasing N, level crossing increase that in Fig.4 clearly shows. 




FIG. 3. Energy levels as a function of the control parameter C in the Hamiltonian (11) for N = 10, 20 bosons wither = 1000, /? = 
-57,5 = 41,7 = 36. 


B. Ground state energy 

The ground state energy is an important observable of phase transition. So, we calculated the ground - state energy 

dE 

,E g . s . and its first derivative , d a c B ' . Fig.5 shows changing of the ground - state energy and its first derivative versus 
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FIG. 4. Same as Fig.3 but with Hamiltonian (12) witha = 1000, /3 = 6.5, 7 = —22. 
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FIG. 5. The ground-state energy (left panel) and derivative of the ground-state energy (right panel) are presented as a function 
of the control parameter C for a system with 10 bosons. 


the control parameter C. Both operators, E g , s , and 


dE a 


dC 


are approximately zero in one phase and different from 


zero in the other phase. Since low N bosons is chosen, it is not possible to distinguish whether the transition is first 
or second order such as done in the even-even case 23,1241. 


C. expectation values of the d-boson number operator 

An appropriate quantal order parameter is: 


{rid) 


W’\iidW 

N 


In order to obtain (rid), we act on the eigenstate, | k; v s vn&LM) 


(rid) 


2 A? - 2C 2 A° + 2fc(y 1 “ 1 - C) 5 
1 -C 2 2 N 


(27) 


Fig.6 shows the expectation values of the d-boson number operator for the lowest states as a function of C control 
parameter for N=10 bosons. Fig.6 displays that the expectation values of the number of d bosons for each J, rid, 
remain approximately constant for C < 0.45 and only begin to change rapidly for C > 0.45. The near constancy of 
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FIG. 6. the expectation values of the d-boson number operator for the lowest states as a function of C control parameter for a 
j = 1/2 particle coupled to a system of (s, d) bosons undergoing a 1/(5) — 0(6) transition (left panel) and single fermion with 
j = 3/2 coupled to a system of ( s,d ) bosons (right panel). 


rid for C < 0.45, is a obvious indication that U (5) dynamical symmetry preserves in this region to a high degree and 
also the rid values change rapidly with C over the range 0.65 < C < 1. It can be seen from Fig .6 that in due to the 
presence of the fermion, the transition is made sharper for states that a j = 1/2 particle coupled to a system of ( s , d) 
bosons undergoing a E/(5) — 0(6) transition (left panel) while is made smoother for a j = 3/2 particle coupled to a 
system of (s,d) bosons (right panel). 


D. energy differences 

Fig.7 displays continues energy differences in terms of control parameter, C, for states that a j = 1/2 and j = 3/2 
particle coupled to a system of (s, d) bosons. Fig.7 shows that during transition from one limit to another exist the 
points that energy is minimum or maximum near the critical point. 


IV. EXPERIMENTAL EVIDENCE 

This section presented the calculated results of low-lying states of the odd -A Ba and Rh isotopes. The results 
include energy levels and the B( E2)values and two-neutron separation energies. 


A. energy spectrum 

Nuclei in the mass regions around A ~ 100 [25], H(| and 130 H3] have transitional characteristics intermediate 
between spherical and gamma-unstable shapes . The theoretical and experimental studies of energy spectra done in 
refs. [2(1 l28l. |20] show Rh and Ba isotopes have 17(5 ) o 0(6) transitional characteristics. .The possible occurrence of 
this symmetry in H 5 Ba has been recently suggested j28| .The negative parity states in the odd - even nuclei Rh are built 
mainly on the 2pi shell model orbit [~26| . The single - particle orbits lgi, 2ds, 2dsand3si establish the positive parity 

2 I | 2 2 2 2 

states in odd-mass Ba isotopes[291]. In this study, a simplifying assumption is made that single particle states are built 
on the j = 1/2 and j = 3/2. We therefore analyze the negative parity states of the odd-proton nuclei, YA~ lm Rh and 
positive parity states of the odd-neutron nuclei, g 6 Ba. In order to obtain energy spectrum and realistic calculation 
for these nuclei, we need to specify Hamiltonian parameters (11), (12). Eigenvalues of these systems are obtained 
by solving Bethe-Ansatz equations with least square fitting processes to experimental data to obtain constants of 
Hamiltonian. The best fits for Hamiltonian’s parameters, namely a , /3 , S and 7 , used in the present work are shown 
in table 1, 2.Tables 3(a,b...e) and 4(a,b...f) are shown calculated energy spectra along with the experimental values 
. Figure .8 and Figure.9 are also shown a comparison between the available experimental levels and the predictions 
of our results for the 1 ° 1 - 109 Rh and ^‘~ 137 Ba isotopes in the low-lying region of spectra. An acceptable degree 
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FIG. 7. Continues energy differences in terms of control parameter, C, for states that a j = 3/2 (top) and j = l/2(bottom) 
particle coupled to a system of (s , d) bosons. 


Table 1. Parameters of Hamiltonian (11) used in the calculation of the 
Rh Isotopes. All parameters are given in keV. 


Nucleus N C a (5 8 7 a 

l^Rh 5 0.0(1 105.99 1.2572 1.(1782 1.2336 194.78 

J fRh 6 0.46 52 -1.053 1.5518 23.217 136.5 

IfRh 7 0.54 70.61 1.6273 8.7337 -0.0331 97.38 

45 7 Rh 8 0.65 196.30 3.5077 -11.734 7.7931 120.48 

H 9 Rh 9 0.7 245.269 2.438 -24.072 -181.02 184.79 


of agreement is obvious between them. We have tried to extract the best set of parameters which reproduce these 
complete spectra with minimum variations. It means that our suggestion to use this transitional Hamiltonian for the 
description of the Rh and Ba isotopic chain would not have any contradiction with other theoretical studies done 
with special hypotheses about mixing of intruder and normal configurations. On the other hand, predictions of our 
model for the control parameter of considered nuclei, C, describe the vibrational, i.e.C = 0, or rotational, namely 
0=1, confirm this mixing of both vibrating and rotating structures in these nuclei when C ~ 0.5 —> 0.65. Fig.10 and 
Fig. 11 display a comparison between the calculated continues energy differences and experimental data for Ba and 
Rh Isotopes, respectively. It can be seen from figs that our results for Ba Isotopes are better than for Rh Isotopes. 


One of the most basic structural predictions of U BF { 5) — O bf ( 6 ) transition is a ) value. The ratio equal to 

2.2 — 2.3 indicates the spectrum of transitional nuclei I HE El- Thus we calculated this quantity for Rh and Ba 
Isotopes. Fig. 12 shows ) prediction values for Rh and Ba Isotopes. For Rh isotopes this value evolves from 

2.1 to 2.8 while Ba isotopes vary of 2.85 to 1.9. Fig.12 displays that ) values for 45 5 Rh and gg 3-13 5 Ba isotopes 

are approximately 2.2 — 2.4. 
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Table 2. Parameters of Hamiltonian (11) and (12) used in the calculation 
of the Ba Isotopes. All parameters are given in keV. 

Nucleus N C a [3 <5 7 a 

ifBa 8 0.78 3.46 -0.0303 -1.1481 19.098 62.48 

IfBa 7 0.8 13.77 -2.685 -8.5776 15.69 94.25 

IfBa 6 0.77 0.578 -3.3839 -3.3319 27.12 128 

IfBa 5 0.68 21.37 -0.0274 16.46 -1.75 119 

IfBa 4 0.65 48.69 3.32 - 36.57 86.84 

IfBa 3 0.75 366 0.194 - 23.047 128.6 


Table 3a. .Energy spectra for Rh isotope 

if Rh r 

A' 


33 exp 

33cal 

(i/2)r 

2 

0 

0 

0 

(3/2)r 

2 

1 

305.5 

267.3 

(s/2)r 

2 

1 

305.5 

273.5 

( 3 / 2)2 

1 

2 

355.3 

665.8 

( 5 / 2)2 

1 

2 

355.3 

671.9 

(7/2)r 

1 

2 

851.4 

704.1 

(9/2)r 

1 

2 

851.4 

715.2 

( 9 / 2)2 

1 

2 

899.3 

821.2 

( 5 / 2 ) 3 - 

2 

1 

996.4 

794.6 

(3/2)1 

1 

2 

1058 

997.9 

(5/2)1 

1 

2 

1058 

971.8 

( 1 / 2)2 

2 

0 

1531 

1472.9 


Table 3b. .Energy spectra for Rh isotope 

if Rh r 

A' 

fd, 

-L^exp 

33cal 

(i/2)r 

3 

0 

0 

0 

(3/2)r 

2 

1 

294.984 

334 

(s/2)r 

2 

2 

357.408 

450.6 

( 1 / 2)2 

3 

0 

803.07 

598.6 

( 3 / 2)2 

2 

2 

803.07 

692.4 

(7/2)r 

2 

2 

847.58 

705.4 

( 5 / 2)2 

2 

2 

880.47 

809 

(9/2)r 

2 

2 

920.1 

915.3 

( 3 / 2)3 

2 

1 

1277.04 

1226.7 

( 13 / 2 ) 3 - 

1 

3 

1637.64 

1576.5 

( 15 / 2 ) 3 - 

1 

4 

2221.2 

2036.3 

( 17 / 2 ) 3 - 

1 

4 

2345.35 

2432.6 

(17/2)J 

0 

5 

2418.6 

2065.5 


Table 3c. .Energy spectra for 

45 5 i?h isotope 

if Rh r 

A' 

Vd 

33exp 

33 cal 

( 1 / 2 ) 3 - 

3 

0 

129.781 

129.8 

(3/2)r 

3 

1 

392.65 

229.2 

( 5 / 2 ) 3 - 

3 

1 

455.61 

511.7 

( 3 / 2)2 

2 

2 

762.11 

811.2 

( 3 / 2)3 

3 

1 

783 

703.6 

( 5 / 2)2 

2 

2 

817 

821.7 

( 7 / 2 ) 3 - 

2 

2 

817 

833.1 

( 5 / 2)3 

2 

3 

866 

868.2 

( 7 / 2)2 

2 

3 

898 

861.5 

( 7 / 2)3 

2 

3 

976 

937.9 

(9/2)r 

2 

2 

976 

1043.5 

( 3 / 2 ) 4 - 

2 

2 

1147 

1312.7 

( 5 / 2)4 

2 

2 

1147 

1311.9 
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Table 3d. .Energy spectra for J ® 7 Rh isotope 


it'Rh r 

A' 

Z4i 

F, 

J-^exp 

E C al 

(i/ 2 )r 

4 

0 

268.36 

268.4 

(3/2)r 

3 

1 

485.66 

345.4 

(s/ 2 )r 

3 

1 

543.84 

410.9 

(9/2)r 

3 

2 

559.97 

423.4 

( 3 / 2)2 

3 

2 

752.55 

818.8 

( 5 / 2)2 

3 

2 

877.75 

870.4 

( 3 / 2)3 

3 

1 

974.44 

921.9 

( 5 / 2)3 

2 

3 

974.44 

958.6 

(7/2)r 

3 

2 

974.44 

780.9 

(3/2)1 

3 

2 

1009.76 

1016.3 

(5/2)1 

3 

2 

1009.76 

1055.3 

(7/2)2 

2 

3 

1251 

1099.5 

(1/2)2 

4 

0 

1334 

1431.2 


Table 3e. .Energy spectra for 

4 g 9 Rh isotope 

i^Rh r 

A 

Vd 

Eexp 

E C al 

(i/2)r 

4 

0 

374.1 

374.1 

(3/2)r 

4 

1 

568.2 

436.6 

(s/2)r 

4 

1 

623.2 

401.5 

( 3 / 2)2 

3 

2 

704.9 

758.1 

( 5 / 2)2 

3 

2 

856.1 

819.6 

( 5 / 2 ) 3 - 

3 

3 

926.9 

827.6 

( 3 / 2 ) 3 - 

4 

1 

1162.3 

1378.8 

(3/2)1 

3 

2 

1214.3 

1003.3 

(5/2)1 

3 

2 

1283.9 

985.6 

(1/2)2 

4 

0 

1631 

1876.6 

( 3 / 2)5 

4 

1 

1631 

1627.4 

( 1 / 2 ) 3 - 

4 

0 

1753 

1703 

(3/2)a 

4 

1 

1753 

1750 


B. B(E 2) transition 


The observables such as electric quadrupole transition probabilities, B(E 2), as well as quadrupole moment ratios 
within the low-lying state provide important information about QPTs. In this section we discuss the calculation of 
E2 transition strengths and compare the results with the available experimental. The electric quadrupole transition 


Table 4a. .Energy spectra for || 7 Ba isotope 


se ''Ba r 

K 


F, 

J-^exp 

E C al 

(l/2tf 

4 

0 

0 

0 

(3/2)+ 

3 

1 

56.1 

51.4 

(5/2)+ 

3 

1 

81 

146.9 

(7/2)+ 

3 

2 

195.1 

265.7 

(3/2)+ 

3 

2 

269.5 

256.3 

(5/2)+ 

3 

2 

269.5 

237.2 

(7/2)+ 

2 

3 

324.1 

376.3 

(7/2)+ 

2 

3 

374.9 

367.1 

(9/2)+ 

3 

2 

415.6 

526.1 

(11/2)+ 

2 

3 

668.9 

723.8 

(11/2)+ 

2 

4 

867.9 

759.4 

(13/2)+ 

2 

3 

963.6 

972.1 

(15/2)+ 

2 

4 

1291.2 

1259.8 

(15/2)+ 

1 

5 

1519.6 

1425.5 

(17/2)+ 

2 

4 

1654.4 

1584.4 
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Table 4b. .Energy spectra for B a isotope 


IfBa r 

K 

Vd 

E 

J-^exp 

E C al 

PT 

3 

0 

0 

0 

(7/2)+ 

2 

2 

8.42 

48.71 

(3/2)+ 

3 

1 

110.56 

185.67 

(3/2)+ 

2 

2 

253.77 

210.76 

(9/2)+ 

2 

2 

263.1 

249.9 

(1/2)4 

3 

0 

278.58 

331.88 

(3/2)+ 

3 

1 

278.58 

228.15 

(5/2)/ 

3 

1 

278.58 

292.83 

(5/2)+ 

2 

2 

318.4 

290.66 

(3/2)+ 

2 

2 

457.03 

483.53 

(3/2)/ 

3 

1 

459.29 

492.44 

(5/2)/ 

2 

3 

459.29 

548.4 

(7/2)/ 

2 

3 

467.3 

589.6 

(H/2)/ 

2 

3 

544.7 

441.257 

(13/2)/ 

2 

3 

864.1 

918.66 


Table 4c. .Energy spectra iorl^Ba isotope 


si VBa r 

K 

Vd 

Eexp 

E C al 

(1/2)/ 

3 

0 

0 

0 

(3/2)/ 

2 

1 

108.08 

70.5 

(3/2)/ 

2 

2 

285.25 

359.2 

(5/2)/ 

2 

1 

316.587 

424.7 

(1/2)/ 

3 

0 

365.16 

487.8 

(3/2)/ 

2 

1 

525.85 

577.3 

(5/2)/ 

2 

2 

525.85 

494.8 

(7/2)/ 

2 

2 

543.11 

638.1 

(3/2)/ 

2 

1 

561.752 

643.5 

(5/2)/ 

1 

3 

561.75 

473.3 

(3/2)/ 

0 

5 

718.78 

573.7 

(5/2)+ 

2 

1 

718.78 

778.6 

(1/2)/ 

1 

3 

719.49 

698.1 

(15/2)/ 

1 

4 

1796.4 

1821.6 

(17/2)/ 

1 

4 

2121.7 

2282.6 


operator f'( E2 '> in odd-A nuclei consists of a bosonic and a fermionic part[ 


With 


rp(E2) _ rji(E 2 ) rji(E 2 ) 


t b E L 2) = g 2 [s+ x d + d + x s]< 2) + x d]L 2) = QbQb , p 


Qb,m = [s+ x J+ d+ x s]J t 2) + ,\'[d + x d/ 2) 


(28) 

(29) 

(30) 


T ( F E2) x ^'] (2) ( 31 ) 

jj' 

Where Qb and Qyy are the boson and fermion quadrupole operator and qB and qf are the effective boson and fermion 
charges[g, [T|.Iii the j = 1/2 case, the E2 transitions are completely determined by the bosonic part of the E2 
operator . The bosonic part have the specific selection rules, where for former term Aud = ±1, \AL\ < 2 and for 
latter Ai/d = 0, ±2, |AL| < 0,4. The reduced electric quadrupole transition rate between the J, —> Jf states is given 
by 1 


B(E 2; cm Ji —> a.fJf) 


|(a / J / ||r^ 2 )||a I J I )| 2 


2 Ji + 1 


(32) 
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Table 4d. .Energy spectra for l^Ba isotope 


si ?Ba r 

K 

I'd 

J-^exp 

E C al 

(1/2 4 

2 

0 

0 

0 

(3/2)+ 

2 

1 

12.3 

64.86 

(5/2)+ 

2 

1 

291.2 

175.93 

(3/2)+ 

1 

2 

302.4 

419.34 

(1/2)4 

1 

2 

539.8 

443.23 

(7/2)4 

1 

2 

577.5 

701.22 

(3/2)+ 

1 

3 

630.6 

705.24 

(5/2)+ 

1 

2 

630.6 

720.84 

(5/2)+ 

1 

3 

676.5 

675.45 

(3/2)+ 

0 

5 

676.5 

726.015 

(5/2)+ 

1 

3 

858.5 

696.45 

(9/2)+ 

1 

2 

883.3 

933.86 

(7/2)+ 

1 

3 

1112.3 

938.98 


Table 4e. .Energy spectra iorl^Ba isotope 


&°Ba J n 

K 

I'd 

Eexp 

E C al 

ww~ 

2 

0 

0 

0 

(1/2 4 

1 

1 

220.954 

200.66 

(5/2)+ 

1 

1 

480.52 

502.31 

(3/2)4 

1 

1 

587.82 

696.93 

(3/2) + 

2 

0 

855 

866.9 

(1/2)4 

1 

2 

910.35 

985.569 

(5/2)4 

1 

2 

979.966 

904.815 

(3/2)+ 

1 

2 

979.969 

1145.37 


The electric quadrupole moment for a state with spin J is given by Q 

= \/¥< 1 ' J| lTtE2> 11 J> (33) 

For evaluating B(E2), We consider eigenstates (14) that the normalization factor obtain as: 


\ T~4 i ^C 2 (^'—+ 4)) , 2(fc— p+^(vd+ j)) 

\ Hp=l 2-~ii=p\ (1-C 2 y k+1 - P )(l-C 2 yi) ' (l-Vk+i-p){^~Vi) ' 




FIG. 8. Comparison between calculated and experimental spectra of positive parity states in Ba Isotopes. The parameters of 
the calculation are given in Tables 2. In the experimental spectra, taken from f2l)l . 




































14 




FIG. 9. Comparison between calculated and experimental spectra of negative parity states in Rh Isotopes. The parameters of 
the calculation are given in Tables 1. In the experimental spectra, taken from l2pl . 




FIG. 10. A comparison between the calculated continues energy differences and experimental data for Ba. In the experimental 
spectra, taken from [21]. 


Unfortunately there are very few experimental data available on electromagnetic properties for odd -mass Ba [29^ and 
Rh [25| isotopes. The values of effective charge ( qB,qf ) are tabulated in Table.5. Table .6 shows experimental and 
calculated values for B(E2) for negative parity states of 45 3 Rh , and positive parity states of 5 Q 5 Ba. The quadrupole 
moments for 45 3 Rh and gg 5 Ba also display in table.7. Table .6 and Table.7 show that in general there are a good 
agreement between the calculated B(E2) values and the quadrupole moments with the experimentally data. The 
values of the control parameter, C, suggest structural changes in nuclear deformation and shape-phase transitions 
in Odd-proton Rh isotopes and Odd-neutron Ba isotopes. Because of the effect of single nucleon on the transition 
and especially critical point, exact selection of the critical point is difficult. With considering this problem, we 
proposed C ~ 0.5 — 0.65 as critical point. So, we conclude from the values of control parameter which has been 
obtained, value and energy differences that 4 ° 5 Rh and H 3 ~ 135 Ba isotopes are as the best candidates for 

U BF ( 5) — O bf { 6 ) transition. 


C. two-neutron separation energies 

Shape phase transitions in nuclei can be studied experimentally by considering the behavior of the ground state 
energies of a series of isotopes, or, more conveniently, the behavior of the two-neutron separation energies, S 2 n 
Q .On the other hand, the ground-state two-neutron separation energies, S 2 m are observables very sensitive to the 
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Table 4f. .Energy spectra ior}^7Ba isotope 


56 B a r 

K 

I'd 

J-sexp 

E C al 

(372)7“ 

2 

0 

0 

0 

(7/2)7 

1 

1 

1252.5 

1093.9 

(5/2)7 

1 

1 

1294 

1124.25 

(3/2)7 

1 

1 

1463.8 

1607 

(1/2)7 

1 

2 

1857 

1782.9 

(1/2)7 

1 

1 

1836.2 

1899.2 

(3/2)+ 

2 

0 

2041 

1940.4 

(5/2)7 

1 

2 

1899 

2063.2 

(3/2)+ 

1 

2 

1899 

1883.5 

(5/2)+ 

1 

2 

2041 

2041.8 

(7/2)7 

1 

2 

2230 

2258 

(9/2)7 

1 

2 

2230 

2248.6 

(7/2)7 

1 

2 

2340 

2263.5 

(7/2)7 

1 

1 

2423.8 

2530 


Table 5. The coefficients of 
T(E 2 )nsed in the present work 
for Rh Rh and l^Ba 
isotopes. 

Nucleus qB(eb) qf(eb) 

45 s Rh 0.461 0 

IfBa 4.7329 -0.7194 


Table 6. B(E2)values for Rh and 5 Q 5 Ba isotopes. The experimental 
data for \^Rh isotope are taken from | 2l l l25l.The experimental data 
Ba isotope are taken from [2l l l31J| . 


for 135 
Ior 56 


Nucleus Ji 
ifRh (3/2)- 

(5/2)7 

(5/2)J 

(5/2)7 

(5/2)7 

(7/2)7 

(9/2)7 

IfBa (1/2)+ 

(5/2)7 

(7/2)7 

(1/2)7 

(3/2)7 

(3/2)7 



Jj 

B(E 2 (e z b- ! ) 

exp. 

calc. 

(i/2)r 

0.109 

0.1172 

(i/2)r 

0.118 

0.1172 

(i/2)r 

0.0044 

0.0044 

(3/2)r 

0.0768 

0.0645 

(s/2)r 

0.015 

0.0097 

(3/2)7 

0.13 

0.1165 

(5/2)7 

0.179 

0.1349 

(3/2)7 

4.6 

3.696 

(1/2)7 

2.65 

2.913 

(3/2)7 

19.9 

14.784 

(3/2)7 

11.7 

14.403 

(3/2)7 

18 

14.785 

(3/2)7 

7 

7.001 


Table 7. Quadrupole moments for Rh and 
isotopes. The experimental data for Rh isotope are 


taken from [21i |2 
isotope are taken from [2. 


.The experimental data for l^Ba 


Nucleus J n 

WHrf, TTF: 


iTRh 

(1/2)7 

0 

0 


(3/2)7 

-0.32 

-0.2588 


(5/2)7 

-0. 41 

- 0.2824 


(7/2)7 

- 

-0.5538 


(9/2)7 

- 

-0.5355 

i 5 Ba 

(3/2)7 

0.146 

0.1509 


(1/2)7 

- 

0.1349 


(5/2)7 

- 

0.5126 
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FIG. 11. A comparison between the calculated continues energy differences and experimental data for Rh. In the experimental 
spectra, taken from [2111 . 




FIG. 12. (E(v d = 2 ))/{E{y d 


1)) prediction values for Rh and Ba Isotopes. In the experimental spectra, taken from (21] 


details of the nuclear structure. The occurrence of continuities in the behavior of two-neutron separation energies 
describe a second-order shape-phase transition between spherical and 7 — unstable rotor limits |l|, 112 .In due to, 
we have investigated the evolution of two-neutron separation energies along the Ba and Rh isotopic chains by both 
experimental and theoretical values, which have been presented in Fig.13. The binding energy as a function of proton 
and neutron number is given by Q 


E B (N n , N v ) = E c + A n N n + A V N„ + ^B„N n (N n - 1) + ^B V N V (N V - 1) + CN n N v + E D (N V , N„) (35) 


Where N W (N„) is the number of proton (neutron) bosons in the valence shell, E c the contribution from the core and 
E B is the contribution to the binding energy due to the deformation. Using Eq. (35), one obtains the following 
relation for the two-neutron separation energyQ: 


S 2n (N n , N v ) = E b (N„, N v ) - E b (N v , N v -l) = A n + BN n + C n N„ + (. E D {N n , N v ) - E D (N„, N v - 1)) (36) 

Using the S 2n empirical values for these isotopic chains (2lJ we have extracted A n + BN n = 15.1857, 20.526 Mev 
and C n = 1.9784, —1.7992 Mev for Ba and Rh, respectively. Then, we obtained the two-neutron separation energies, 
which are shown in Fig.13, together with the experimental values. It can be seen from Fig.13 that exist continuities 
(linear variation) in the behavior of two-neutron separation energies thus the phase transition for Ba and Rh isotopic 
chains is of second order .Our result Confirmed the predictions of done in refs. sm, where they suggest a linear 
variation of S 2n with respect to the neutron number for the 17(5) — SO(6) transitional region. 
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FIG. 13. A comparison between theoretical and experimental two neutron separation energies, S 2 n (in keV) for Ba isotopes 
(left panel) and Rh isotopes (right panel). Experimental data from [2l| . 


V. CONCLUSIONS 


In this paper, we have analyzed transition from spherical to 7 — unstable shapes in odd - A nuclei. Key observables of 
phase transition such as level crossing, ground-state energy and derivative of the ground-state energy and expectation 
values of the d-boson number operator have calculated. We have presented experimental evidence for the U (5) — 0(6) 
transition for negative parity states of the 45 - Rh isotopic chain and positive parity states of 56 ' - Ba isotopic 
chain, and performed an analysis for these isotopes via a SU(l,l)-based Hamiltonian. The results indicate that the 
energy spectra of the Rh and Ba isotopes can be reproduced quite well. The calculated B(E2)values and two neutron 
separation energies are agreements with the available experimental data. Our results show that Rh isotopes have 
gamma-unstable rotor features but the vibrational character is dominant while a dominancy of dynamical symmetry 
0(6) exist for Ba isotopic chain and also 4 g 5 Rh and b Ba isotopes are as the best candidates for U BF (5) — O bf (6) 
transition. 
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